Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data. IMPORTANCE: Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Targeting TRIP13 in favorable histology Wilms tumor with nuclear export inhibitors synergizes with doxorubicin.

Wilms tumor (WT) is the most common renal malignancy of childhood. Despite improvements in the overall survival, relapse occurs in ~15% of patients with favorable histology WT (FHWT). Half of these patients will succumb to their disease. Identifying novel targeted therapies remains challenging in part due to the lack of faithful preclinical in vitro models. Here we establish twelve patient-derived WT cell lines and demonstrate that these models faithfully recapitulate WT biology using genomic and transcriptomic techniques. We then perform loss-of-function screens to identify the nuclear export gene, XPO1, as a vulnerability. We find that the FDA approved XPO1 inhibitor, KPT-330, suppresses TRIP13 expression, which is required for survival. We further identify synergy between KPT-330 and doxorubicin, a chemotherapy used in high-risk FHWT. Taken together, we identify XPO1 inhibition with KPT-330 as a potential therapeutic option to treat FHWTs and in combination with doxorubicin, leads to durable remissions in vivo.

Single cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10x Genomics Chromium single cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human moDCs infected with ZIKV at the single cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN dependent and independent genes (antiviral module). We modeled the ZIKV specific antiviral state at the protein level leveraging experimentally derived protein-interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per cell basis with experimental protein interaction data. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool to gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Integrated single-nuclei and spatial transcriptomic analysis reveals propagation of early acute vein harvest and distension injury signaling pathways following arterial implantation.

Background: Vein graft failure (VGF) following cardiovascular bypass surgery results in significant patient morbidity and cost to the healthcare system. Vein graft injury can occur during autogenous vein harvest and preparation, as well as after implantation into the arterial system, leading to the development of intimal hyperplasia, vein graft stenosis, and, ultimately, bypass graft failure. While previous studies have identified maladaptive pathways that occur shortly after implantation, the specific signaling pathways that occur during vein graft preparation are not well defined and may result in a cumulative impact on VGF. We, therefore, aimed to elucidate the response of the vein conduit wall during harvest and following implantation, probing the key maladaptive pathways driving graft failure with the overarching goal of identifying therapeutic targets for biologic intervention to minimize these natural responses to surgical vein graft injury. Methods: Employing a novel approach to investigating vascular pathologies, we harnessed both single-nuclei RNA-sequencing (snRNA-seq) and spatial transcriptomics (ST) analyses to profile the genomic effects of vein grafts after harvest and distension, then compared these findings to vein grafts obtained 24 hours after carotid-cartoid vein bypass implantation in a canine model (n=4). Results: Spatial transcriptomic analysis of canine cephalic vein after initial conduit harvest and distention revealed significant enrichment of pathways ( P < 0.05) involved in the activation of endothelial cells (ECs), fibroblasts (FBs), and vascular smooth muscle cells (VSMCs), namely pathways responsible for cellular proliferation and migration and platelet activation across the intimal and medial layers, cytokine signaling within the adventitial layer, and extracellular matrix (ECM) remodeling throughout the vein wall. Subsequent snRNA-seq analysis supported these findings and further unveiled distinct EC and FB subpopulations with significant upregulation ( P < 0.00001) of markers related to endothelial injury response and cellular activation of ECs, FBs, and VSMCs. Similarly, in vein grafts obtained 24 hours after arterial bypass, there was an increase in myeloid cell, protomyofibroblast, injury-response EC, and mesenchymal-transitioning EC subpopulations with a concomitant decrease in homeostatic ECs and fibroblasts. Among these markers were genes previously implicated in vein graft injury, including VCAN (versican), FBN1 (fibrillin-1), and VEGFC (vascular endothelial growth factor C), in addition to novel genes of interest such as GLIS3 (GLIS family zinc finger 3) and EPHA3 (ephrin-A3). These genes were further noted to be driving the expression of genes implicated in vascular remodeling and graft failure, such as IL-6 , TGFBR1 , SMAD4 , and ADAMTS9. By integrating the ST and snRNA-seq datasets, we highlighted the spatial architecture of the vein graft following distension, wherein activated and mesenchymal-transitioning ECs, myeloid cells, and FBs were notably enriched in the intima and media of distended veins. Lastly, intercellular communication network analysis unveiled the critical roles of activated ECs, mesenchymal transitioning ECs, protomyofibroblasts, and VSMCs in upregulating signaling pathways associated with cellular proliferation (MDK, PDGF, VEGF), transdifferentiation (Notch), migration (ephrin, semaphorin), ECM remodeling (collagen, laminin, fibronectin), and inflammation (thrombospondin), following distension. Conclusions: Vein conduit harvest and distension elicit a prompt genomic response facilitated by distinct cellular subpopulations heterogeneously distributed throughout the vein wall. This response was found to be further exacerbated following vein graft implantation, resulting in a cascade of maladaptive gene regulatory networks. Together, these results suggest that distension initiates the upregulation of pathological pathways that may ultimately contribute to bypass graft failure and presents potential early targets warranting investigation for targeted therapies. This work highlights the first applications of single-nuclei and spatial transcriptomic analyses to investigate venous pathologies, underscoring the utility of these methodologies and providing a foundation for future investigations.

Indoximod-based chemo-immunotherapy for pediatric brain tumors: A first-in-children phase I trial.

BACKGROUND: Recurrent brain tumors are the leading cause of cancer death in children. Indoleamine 2,3-dioxygenase (IDO) is a targetable metabolic checkpoint that, in preclinical models, inhibits anti-tumor immunity following chemotherapy. METHODS: We conducted a phase I trial (NCT02502708) of the oral IDO-pathway inhibitor indoximod in children with recurrent brain tumors or newly diagnosed diffuse intrinsic pontine glioma (DIPG). Separate dose-finding arms were performed for indoximod in combination with oral temozolomide (200 mg/m2/day x 5 days in 28-day cycles), or with palliative conformal radiation. Blood samples were collected at baseline and monthly for single-cell RNA-sequencing with paired single-cell T cell receptor sequencing. RESULTS: Eighty-one patients were treated with indoximod-based combination therapy. Median follow-up was 52 months (range 39-77 months). Maximum tolerated dose was not reached, and the pediatric dose of indoximod was determined as 19.2 mg/kg/dose, twice daily. Median overall survival was 13.3 months (n = 68, range 0.2-62.7) for all patients with recurrent disease and 14.4 months (n = 13, range 4.7-29.7) for DIPG. The subset of n = 26 patients who showed evidence of objective response (even a partial or mixed response) had over 3-fold longer median OS (25.2 months, range 5.4-61.9, p = 0.006) compared to n = 37 nonresponders (7.3 months, range 0.2-62.7). Four patients remain free of active disease longer than 36 months. Single-cell sequencing confirmed emergence of new circulating CD8 T cell clonotypes with late effector phenotype. CONCLUSIONS: Indoximod was well tolerated and could be safely combined with chemotherapy and radiation. Encouraging preliminary evidence of efficacy supports advancing to Phase II/III trials for pediatric brain tumors.

The Simple prEservatioN of Single cElls method for cryopreservation enables the generation of single-cell immune profiles from whole blood.

Introduction: Current multistep methods utilized for preparing and cryopreserving single-cell suspensions from blood samples for single-cell RNA sequencing (scRNA-seq) are time-consuming, requiring trained personnel and special equipment, so limiting their clinical adoption. We developed a method, Simple prEservatioN of Single cElls (SENSE), for single-step cryopreservation of whole blood (WB) along with granulocyte depletion during single-cell assay, to generate high quality single-cell profiles (SCP). Methods: WB was cryopreserved using the SENSE method and peripheral blood mononuclear cells (PBMCs) were isolated and cryopreserved using the traditional density-gradient method (PBMC method) from the same blood sample (n=6). The SCPs obtained from both methods were processed using a similar pipeline and quality control parameters. Further, entropy calculation, differential gene expression, and cellular communication analysis were performed to compare cell types and subtypes from both methods. Results: Highly viable (86.3 ± 1.51%) single-cell suspensions (22,353 cells) were obtained from the six WB samples cryopreserved using the SENSE method. In-depth characterization of the scRNA-seq datasets from the samples processed with the SENSE method yielded high-quality profiles of lymphoid and myeloid cell types which were in concordance with the profiles obtained with classical multistep PBMC method processed samples. Additionally, the SENSE method cryopreserved samples exhibited significantly higher T-cell enrichment, enabling deeper characterization of T-cell subtypes. Overall, the SENSE and PBMC methods processed samples exhibited transcriptional, and cellular communication network level similarities across cell types with no batch effect except in myeloid lineage cells. Discussion: Comparative analysis of scRNA-seq datasets obtained with the two cryopreservation methods i.e., SENSE and PBMC methods, yielded similar cellular and molecular profiles, confirming the suitability of the former method's incorporation in clinics/labs for cryopreserving and obtaining high-quality single-cells for conducting critical translational research.

Single-cell RNA sequencing distinctly characterizes the wide heterogeneity in pediatric mixed phenotype acute leukemia.

BACKGROUND: Mixed phenotype acute leukemia (MPAL), a rare subgroup of leukemia characterized by blast cells with myeloid and lymphoid lineage features, is difficult to diagnose and treat. A better characterization of MPAL is essential to understand the subtype heterogeneity and how it compares with acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). Therefore, we performed single-cell RNA sequencing (scRNAseq) on pediatric MPAL bone marrow (BM) samples to develop a granular map of the MPAL blasts and microenvironment landscape. METHODS: We analyzed over 40,000 cells from nine pediatric MPAL BM samples to generate a single-cell transcriptomic landscape of B/myeloid (B/My) and T/myeloid (T/My) MPAL. Cells were clustered using unsupervised single-cell methods, and malignant blast and immune clusters were annotated. Differential expression analysis was performed to identify B/My and T/My MPAL blast-specific signatures by comparing transcriptome profiles of MPAL with normal BM, AML, and ALL. Gene set enrichment analysis (GSEA) was performed, and significantly enriched pathways were compared in MPAL subtypes. RESULTS: B/My and T/My MPAL blasts displayed distinct blast signatures. Transcriptomic analysis revealed that B/My MPAL profile overlaps with B-ALL and AML samples. Similarly, T/My MPAL exhibited overlap with T-ALL and AML samples. Genes overexpressed in both MPAL subtypes' blast cells compared to AML, ALL, and healthy BM included MAP2K2 and CD81. Subtype-specific genes included HBEGF for B/My and PTEN for T/My. These marker sets segregated bulk RNA-seq AML, ALL, and MPAL samples based on expression profiles. Analysis comparing T/My MPAL to ETP, near-ETP, and non-ETP T-ALL, showed that T/My MPAL had greater overlap with ETP-ALL cases. Comparisons among MPAL subtypes between adult and pediatric samples showed analogous transcriptomic landscapes of corresponding subtypes. Transcriptomic differences were observed in the MPAL samples based on response to induction chemotherapy, including selective upregulation of the IL-16 pathway in relapsed samples. CONCLUSIONS: We have for the first time described the single-cell transcriptomic landscape of pediatric MPAL and demonstrated that B/My and T/My MPAL have distinct scRNAseq profiles from each other, AML, and ALL. Differences in transcriptomic profiles were seen based on response to therapy, but larger studies will be needed to validate these findings.

Single-cell analysis reveals altered tumor microenvironments of relapse- and remission-associated pediatric acute myeloid leukemia.

Acute myeloid leukemia (AML) microenvironment exhibits cellular and molecular differences among various subtypes. Here, we utilize single-cell RNA sequencing (scRNA-seq) to analyze pediatric AML bone marrow (BM) samples from diagnosis (Dx), end of induction (EOI), and relapse timepoints. Analysis of Dx, EOI scRNA-seq, and TARGET AML RNA-seq datasets reveals an AML blasts-associated 7-gene signature (CLEC11A, PRAME, AZU1, NREP, ARMH1, C1QBP, TRH), which we validate on independent datasets. The analysis reveals distinct clusters of Dx relapse- and continuous complete remission (CCR)-associated AML-blasts with differential expression of genes associated with survival. At Dx, relapse-associated samples have more exhausted T cells while CCR-associated samples have more inflammatory M1 macrophages. Post-therapy EOI residual blasts overexpress fatty acid oxidation, tumor growth, and stemness genes. Also, a post-therapy T-cell cluster associated with relapse samples exhibits downregulation of MHC Class I and T-cell regulatory genes. Altogether, this study deeply characterizes pediatric AML relapse- and CCR-associated samples to provide insights into the BM microenvironment landscape.

Taxane chemotherapy induces stromal injury that leads to breast cancer dormancy escape.

A major cause of cancer recurrence following chemotherapy is cancer dormancy escape. Taxane-based chemotherapy is standard of care in breast cancer treatment aimed at killing proliferating cancer cells. Here, we demonstrate that docetaxel injures stromal cells, which release protumor cytokines, IL-6 and granulocyte colony stimulating factor (G-CSF), that in turn invoke dormant cancer outgrowth both in vitro and in vivo. Single-cell transcriptomics shows a reprogramming of awakened cancer cells including several survival cues such as stemness, chemoresistance in a tumor stromal organoid (TSO) model, as well as an altered tumor microenvironment (TME) with augmented protumor immune signaling in a syngeneic mouse breast cancer model. IL-6 plays a role in cancer cell proliferation, whereas G-CSF mediates tumor immunosuppression. Pathways and differential expression analyses confirmed MEK as the key regulatory molecule in cancer cell outgrowth and survival. Antibody targeting of protumor cytokines (IL-6, G-CSF) or inhibition of cytokine signaling via MEK/ERK pathway using selumetinib prior to docetaxel treatment prevented cancer dormancy outgrowth suggesting a novel therapeutic strategy to prevent cancer recurrence.

Pediatric T-cell acute lymphoblastic leukemia blast signature and MRD associated immune environment changes defined by single cell transcriptomics analysis.

Different driver mutations and/or chromosomal aberrations and dysregulated signaling interactions between leukemia cells and the immune microenvironment have been implicated in the development of T-cell acute lymphoblastic leukemia (T-ALL). To better understand changes in the bone marrow microenvironment and signaling pathways in pediatric T-ALL, bone marrows collected at diagnosis (Dx) and end of induction therapy (EOI) from 11 patients at a single center were profiled by single cell transcriptomics (10 Dx, 5 paired EOI, 1 relapse). T-ALL blasts were identified by comparison with healthy bone marrow cells. T-ALL blast-associated gene signature included SOX4, STMN1, JUN, HES4, CDK6, ARMH1 among the most significantly overexpressed genes, some of which are associated with poor prognosis in children with T-ALL. Transcriptome profiles of the blast cells exhibited significant inter-patient heterogeneity. Post induction therapy expression profiles of the immune cells revealed significant changes. Residual blast cells in MRD+ EOI samples exhibited significant upregulation (P < 0.01) of PD-1 and RhoGDI signaling pathways. Differences in cellular communication were noted in the presence of residual disease in T cell and hematopoietic stem cell compartments in the bone marrow. Together, these studies generate new insights and expand our understanding of the bone marrow landscape in pediatric T-ALL.

DCUN1D1 Is an Essential Regulator of Prostate Cancer Proliferation and Tumour Growth That Acts through Neddylation of Cullin 1, 3, 4A and 5 and Deregulation of Wnt/Catenin Pathway.

Defective in cullin neddylation 1 domain containing 1 (DCUN1D1) is an E3 ligase for the neddylation, a post-translational process similar to and occurring in parallel to ubiquitin proteasome pathway. Although established as an oncogene in a variety of squamous cell carcinomas, the precise role of DCUN1D1 in prostate cancer (PCa) has not been previously explored thoroughly. Here, we investigated the role of DCUN1D1 in PCa and demonstrated that DCUN1D1 is upregulated in cell lines as well as human tissue samples. Inhibition of DCUN1D1 significantly reduced PCa cell proliferation and migration and remarkably inhibited xenograft formation in mice. Applying both genomics and proteomics approaches, we provide novel information about the DCUN1D1 mechanism of action. We identified CUL3, CUL4B, RBX1, CAND1 and RPS19 proteins as DCUN1D1 binding partners. Our analysis also revealed the dysregulation of genes associated with cellular growth and proliferation, developmental, cell death and cancer pathways and the WNT/ß-catenin pathway as potential mechanisms. Inhibition of DCUN1D1 leads to the inactivation of ß-catenin through its phosphorylation and degradation which inhibits the downstream action of ß-catenin, reducing its interaction with Lef1 in the Lef1/TCF complex that regulates Wnt target gene expression. Together our data point to an essential role of the DCUN1D1 protein in PCa which can be explored for potential targeted therapy.

The 'omics of obesity in B-cell acute lymphoblastic leukemia.

The obesity pandemic currently affects more than 70 million Americans and more than 650 million individuals worldwide. In addition to increasing susceptibility to pathogenic infections (eg, SARS-CoV-2), obesity promotes the development of many cancer subtypes and increases mortality rates in most cases. We and others have demonstrated that, in the context of B-cell acute lymphoblastic leukemia (B-ALL), adipocytes promote multidrug chemoresistance. Furthermore, others have demonstrated that B-ALL cells exposed to the adipocyte secretome alter their metabolic states to circumvent chemotherapy-mediated cytotoxicity. To better understand how adipocytes impact the function of human B-ALL cells, we used a multi-omic RNA-sequencing (single-cell and bulk transcriptomic) and mass spectroscopy (metabolomic and proteomic) approaches to define adipocyte-induced changes in normal and malignant B cells. These analyses revealed that the adipocyte secretome directly modulates programs in human B-ALL cells associated with metabolism, protection from oxidative stress, increased survival, B-cell development, and drivers of chemoresistance. Single-cell RNA sequencing analysis of mice on low- and high-fat diets revealed that obesity suppresses an immunologically active B-cell subpopulation and that the loss of this transcriptomic signature in patients with B-ALL is associated with poor survival outcomes. Analyses of sera and plasma samples from healthy donors and those with B-ALL revealed that obesity is associated with higher circulating levels of immunoglobulin-associated proteins, which support observations in obese mice of altered immunological homeostasis. In all, our multi-omics approach increases our understanding of pathways that may promote chemoresistance in human B-ALL and highlight a novel B-cell-specific signature in patients associated with survival outcomes.

Cross center single-cell RNA sequencing study of the immune microenvironment in rapid progressing multiple myeloma.

Despite advancements in understanding the pathophysiology of Multiple Myeloma (MM), the cause of rapid progressing disease in a subset of patients is still unclear. MM's progression is facilitated by complex interactions with the surrounding bone marrow (BM) cells, forming a microenvironment that supports tumor growth and drug resistance. Understanding the immune microenvironment is key to identifying factors that promote rapid progression of MM. To accomplish this, we performed a multi-center single-cell RNA sequencing (scRNA-seq) study on 102,207 cells from 48 CD138- BM samples collected at the time of disease diagnosis from 18 patients with either rapid progressing (progression-free survival (PFS) < 18 months) or non-progressing (PFS > 4 years) disease. Comparative analysis of data from three centers demonstrated similar transcriptome profiles and cell type distributions, indicating subtle technical variation in scRNA-seq, opening avenues for an expanded multicenter trial. Rapid progressors depicted significantly higher enrichment of GZMK+ and TIGIT+ exhausted CD8+ T-cells (P = 0.022) along with decreased expression of cytolytic markers (PRF1, GZMB, GNLY). We also observed a significantly higher enrichment of M2 tolerogenic macrophages in rapid progressors and activation of pro-proliferative signaling pathways, such as BAFF, CCL, and IL16. On the other hand, non-progressive patients depicted higher enrichment for immature B Cells (i.e., Pre/Pro B cells), with elevated expression for markers of B cell development (IGLL1, SOX4, DNTT). This multi-center study identifies the enrichment of various pro-tumorigenic cell populations and pathways in those with rapid progressing disease and further validates the robustness of scRNA-seq data generated at different study centers.

Therapeutic implications of mitochondrial stress-induced proteasome inhibitor resistance in multiple myeloma.

The connections between metabolic state and therapy resistance in multiple myeloma (MM) are poorly understood. We previously reported that electron transport chain (ETC) suppression promotes sensitivity to the BCL-2 antagonist venetoclax. Here, we show that ETC suppression promotes resistance to proteasome inhibitors (PIs). Interrogation of ETC-suppressed MM reveals integrated stress response-dependent suppression of protein translation and ubiquitination, leading to PI resistance. ETC and protein translation gene expression signatures from the CoMMpass trial are down-regulated in patients with poor outcome and relapse, corroborating our in vitro findings. ETC-suppressed MM exhibits up-regulation of the cystine-glutamate antiporter SLC7A11, and analysis of patient single-cell RNA-seq shows that clusters with low ETC gene expression correlate with higher SLC7A11 expression. Furthermore, erastin or venetoclax treatment diminishes mitochondrial stress-induced PI resistance. In sum, our work demonstrates that mitochondrial stress promotes PI resistance and underscores the need for implementing combinatorial regimens in MM cognizant of mitochondrial metabolic state.

Transcriptomics and systems biology identify non-antibiotic drugs for the treatment of ocular bacterial infection.

Increasing antibiotic resistance among ocular pathogens often results in treatment failure for blinding infections such as endophthalmitis. Hence, newer therapeutics is needed to combat multidrug-resistant infections. Here, we show a drug repurposing approach using a connectivity map based on temporal transcriptomics of Staphylococcus aureus (SA) infected mouse retina. The analysis predicted three non-antibiotic drugs, Dequalinium chloride (DC), Clofilium tosylate (CT), and Glybenclamide (Glb) which reversed the SA infection signatures. Predicted drugs exhibited anti-inflammatory properties in human retinal cells against sensitive and resistant strains of SA. Intravitreal administration of all drugs reduced intraocular inflammation in SA-infected mouse eyes while DC and CT also reduced bacterial burden. Drug treatment improved visual function coinciding with reduced Caspase-3 mediated retinal cell death. Importantly, all drugs exhibited synergy with vancomycin in improving disease outcomes. Overall, our study identified three non-antibiotic drugs and demonstrated their therapeutic and prophylactic efficacies in ameliorating intraocular bacterial infection.

Cellular and immunometabolic mechanisms of inflammation in depression: Preliminary findings from single cell RNA sequencing and a tribute to Bruce McEwen.

Inflammation is associated with symptoms of anhedonia, a core feature of major depression (MD). We have shown that MD patients with high inflammation as measured by plasma C-reactive protein (CRP) and anhedonia display gene signatures of metabolic reprograming (e.g., shift to glycolysis) necessary to sustain cellular immune activation. To gain preliminary insight into the immune cell subsets and transcriptomic signatures that underlie increased inflammation and its relationship with behavior in MD at the single-cell (sc) level, herein we conducted scRNA-Seq on peripheral blood mononuclear cells from a subset of medically-stable, unmedicated MD outpatients. Three MD patients with high CRP (>3 mg/L) before and two weeks after anti-inflammatory challenge with the tumor necrosis factor antagonist infliximab and three patients with low CRP (≤3 mg/L) were studied. Cell clusters were identified using a Single Cell Wizard pipeline, followed by pathway analysis. CD14+ and CD16+ monocytes were more abundant in MD patients with high CRP and were reduced by 29% and 55% respectively after infliximab treatment. Within CD14+ and CD16+ monocytes, genes upregulated in high CRP patients were enriched for inflammatory (phagocytosis, complement, leukocyte migration) and immunometabolic (hypoxia-inducible factor [HIF]-1, aerobic glycolysis) pathways. Shifts in CD4+ T cell subsets included ∼30% and ∼10% lower abundance of CD4+ central memory (TCM) and naïve cells and ∼50% increase in effector memory-like (TEM-like) cells in high versus low CRP patients. TCM cells of high CRP patients displayed downregulation of the oxidative phosphorylation (OXPHOS) pathway, a main energy source in this cell type. Following infliximab, changes in the number of CD14+ monocytes and CD4+ TEM-like cells predicted improvements in anhedonia scores (r = 1.0, p < 0.001). In sum, monocytes and CD4+ T cells from MD patients with increased inflammation exhibited immunometabolic reprograming in association with symptoms of anhedonia. These findings are the first step toward determining the cellular and molecular immune pathways associated with inflammatory phenotypes in MD, which may lead to novel immunomodulatory treatments of psychiatric illnesses with increased inflammation.

Baricitinib attenuates the proinflammatory phase of COVID-19 driven by lung-infiltrating monocytes.

SARS-CoV-2-infected subjects are generally asymptomatic during initial viral replication but may suffer severe immunopathology after the virus has receded and monocytes have infiltrated the airways. In bronchoalveolar lavage fluid from severe COVID-19 patients, monocytes express mRNA encoding inflammatory mediators and contain SARS-CoV-2 transcripts. We leverage a human small airway model of infection and inflammation, whereby primary blood monocytes transmigrate across SARS-CoV-2-infected lung epithelium to characterize viral burden, gene expression, and inflammatory mediator secretion by epithelial cells and monocytes. In this model, lung-infiltrating monocytes acquire SARS-CoV-2 from the epithelium and upregulate expression and secretion of inflammatory mediators, mirroring in vivo data. Combined use of baricitinib (Janus kinase inhibitor) and remdesivir (nucleoside analog) enhances antiviral signaling and viral clearance by SARS-CoV-2-positive monocytes while decreasing secretion of proneutrophilic mediators associated with acute respiratory distress syndrome. These findings highlight the role of lung-infiltrating monocytes in COVID-19 pathogenesis and their importance as a therapeutic target.

Comprehensive Genomic Profiling of High-Risk Pediatric Cancer Patients Has a Measurable Impact on Clinical Care.

PURPOSE: Profiling of pediatric cancers through deep sequencing of large gene panels and whole exomes is rapidly being adopted in many clinical settings. However, the most impactful approach to genomic profiling of pediatric cancers remains to be defined. METHODS: We conducted a prospective precision medicine trial, using whole-exome sequencing of tumor and germline tissue and whole-transcriptome sequencing (RNA Seq) of tumor tissue to characterize the mutational landscape of 127 tumors from 126 unique patients across the spectrum of pediatric brain tumors, hematologic malignancies, and extracranial solid tumors. RESULTS: We identified somatic tumor alterations in 121/127 (95.3%) tumor samples and identified cancer predisposition syndromes on the basis of known pathogenic or likely pathogenic germline mutations in cancer predisposition genes in 9/126 patients (7.1%). Additionally, we developed a novel scoring system for measuring the impact of tumor and germline sequencing, encompassing therapeutically relevant genomic alterations, cancer-related germline findings, recommendations for treatment, and refinement of risk stratification or prognosis. At least one impactful finding from the genomic results was identified in 108/127 (85%) samples sequenced. A recommendation to consider a targeted agent was provided for 82/126 (65.1%) patients. Twenty patients ultimately received therapy with a molecularly targeted agent, representing 24% of those who received a targeted agent recommendation and 16% of the total cohort. CONCLUSION: Paired tumor/normal whole-exome sequencing and tumor RNA Seq of de novo or relapsed/refractory tumors was feasible and clinically impactful in high-risk pediatric cancer patients.

Autocrine Canonical Wnt Signaling Primes Noncanonical Signaling through ROR1 in Metastatic Castration-Resistant Prostate Cancer.

Wnt signaling driven by genomic alterations in genes including APC and CTNNB, which encodes ß-catenin, have been implicated in prostate cancer development and progression to metastatic castration-resistant prostate cancer (mCRPC). However, nongenomic drivers and downstream effectors of Wnt signaling in prostate cancer and the therapeutic potential of targeting this pathway in prostate cancer have not been fully established. Here we analyzed Wnt/ß-catenin signaling in prostate cancer and identified effectors distinct from those found in other tissues, including aryl hydrocarbon receptor and RUNX1, which are linked to stem cell maintenance, and ROR1, a noncanonical Wnt5a coreceptor. Wnt/ß-catenin signaling-mediated increases in ROR1 enhanced noncanonical responses to Wnt5a. Regarding upstream drivers, APC genomic loss, but not its epigenetic downregulation commonly observed in prostate cancer, was strongly associated with Wnt/ß-catenin pathway activation in clinical samples. Tumor cell upregulation of the Wnt transporter Wntless (WLS) was strongly associated with Wnt/ß-catenin pathway activity in primary prostate cancer but also associated with both canonical and noncanonical Wnt signaling in mCRPC. IHC confirmed tumor cell WLS expression in primary prostate cancer and mCRPC, and patient-derived prostate cancer xenografts expressing WLS were responsive to treatment with Wnt synthesis inhibitor ETC-1922159. These findings reveal that Wnt/ß-catenin signaling in prostate cancer drives stem cell maintenance and invasion and primes for noncanonical Wnt signaling through ROR1. They further show that autocrine Wnt production is a nongenomic driver of canonical and noncanonical Wnt signaling in prostate cancer, which can be targeted with Wnt synthesis inhibitors to suppress tumor growth. SIGNIFICANCE: This work provides fundamental insights into Wnt signaling and prostate cancer cell biology and indicates that a subset of prostate cancer driven by autocrine Wnt signaling is sensitive to Wnt synthesis inhibitors.

Survival Genie, a web platform for survival analysis across pediatric and adult cancers.

The genomics data-driven identification of gene signatures and pathways has been routinely explored for predicting cancer survival and making decisions related to targeted treatments. A large number of packages and tools have been developed to correlate gene expression/mutations to the clinical outcome but lack the ability to perform such analysis based on pathways, gene sets, and gene ratios. Furthermore, in this single-cell omics era, the cluster markers from cancer single-cell transcriptomics studies remain an underutilized prognostic option. Additionally, no bioinformatics online tool evaluates the associations between the enrichment of canonical cell types and survival across cancers. Here we have developed Survival Genie, a web tool to perform survival analysis on single-cell RNA-seq (scRNA-seq) data and a variety of other molecular inputs such as gene sets, genes ratio, tumor-infiltrating immune cells proportion, gene expression profile scores, and tumor mutation burden. For a comprehensive analysis, Survival Genie contains 53 datasets of 27 distinct malignancies from 11 different cancer programs related to adult and pediatric cancers. Users can upload scRNA-seq data or gene sets and select a gene expression partitioning method (i.e., mean, median, quartile, cutp) to determine the effect of expression levels on survival outcomes. The tool provides comprehensive results including box plots of low and high-risk groups, Kaplan-Meier plots with univariate Cox proportional hazards model, and correlation of immune cell enrichment and molecular profile. The analytical options and comprehensive collection of cancer datasets make Survival Genie a unique resource to correlate gene sets, pathways, cellular enrichment, and single-cell signatures to clinical outcomes to assist in developing next-generation prognostic and therapeutic biomarkers. Survival Genie is open-source and available online at https://bbisr.shinyapps.winship.emory.edu/SurvivalGenie/ .